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CN ■ Abstract 

A finite-size scaling theory for the cpf model is derived using renor- 
malization group methods. Particular attention is paid to the partition 
O ; function zeroes, in terms of which all thermodynamic observables can 

O . be expressed. While the leading scaling behaviour is identical to that 

of mean field theory, there exist multiplicative logarithmic corrections 
too. A non-perturbative test of these formulae in the form of a high 
^ . precision Monte Carlo analysis reveals good quantitative agreement 

H I with the analytical predictions. 
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1 Introduction 

Above one dimension, lattice 0^ theory is known to possess a second order 
phase transition separating an ordered phase from a disordered one. The con- 
tinuum parameterization of the field theory is defined at this phase transition. 
There exist rigorous proofs regarding the trivial (Gaussian) and interactive 
nature of the continuum theory in d > 4 and d < 4 dimensions respectively 
[H]. Although a rigorous proof is still lacking, it is believed that in d = 4 
the theory is also trivial. Nonetheless, the theory may be useful as one with 
an effective interaction — valid below some momentum cutoff A. In this 
context, triviality of the theory means that the (leading) critical exponents 
at the phase transition are identical to those of the Gaussian model which 
describes free bosons. The renormalization group (RG) approach predicts 
logarithmic violations of the mean field scaling relations in four dimensions 
[|, H . It has been stated @] that the existence of such logarithmic corrections 
to any mean field scaling relation implies triviality. Hence the importance 
of the study of these corrections and the primary motivation for the present 
work. 

The layout of this paper is as follows: 

In sect. 2 the perturbative RG is applied to the single component 0^ theory, 
and a finite-size scaling (FSS) theory is developed with particular emphasis on 
four dimensions. The approach to criticality from within both the symmetric 
and broken phases is studied. 

The lattice version of the model is then discussed in sect. 3. The con- 
cept of partition function zeroes as an alternative way to view the onset of 
a phase transition is recalled. Formulae describing the FSS behaviour of 
these zeroes are derived. These are then used to derive the FSS formulae of 
thermodynamic functions. 

The RG and FSS equations predict logarithmic violations of the mean 
field theory in four dimensions. These solutions are based on perturbation 
theory and have to be tested in an independent approach. 

In sect. 4 details of our numerical simulations are discussed and the results 
regarding logarithmic scaling corrections presented. These are consistent 
with our RG predictions. We conclude in sect. 5. 

We should mention here that some of the background material is not 
new. It has been included to keep the presentation as self contained as 
possible. Our treatment of the FSS behaviour of the Lee-Yang zeroes (as 



well as the consequent FSS behaviour of the thermodynamic functions) in 
four dimensions is, however, new. Some of our numerical results on Fisher 
zeroes have been presented earlier |]^. 

2 Renormalization group and finite-size scal- 
ing 

In the single component version of (f)^ theory with quadratic composite fields, 
the Hamiltonian density in rf-dimensional Euclidean space-time continuum 

may be written as p 

n = \{V<pr + ^<P\x) + ^^<p\x) - H{x)<p{x) - '-^<P\x). (2.1) 

Here mo is the bare mass of the bosons described by the theory, and go is 
the bare interaction coupling. H{x) and t{x) are the sources for the fields (f) 
and the composite (j)"^ fields respectively. The generating functionals Z[H,t] 
and W [H, t] are defined by 

Z[H,t\ = e^[^'*] = C fl[d(t)ix)e~I^'''^, (2.2) 

and the constant C is chosen such that 

H^[0,0]=0. (2.3) 

The function conjugate to H{x) is 

M{x,t) = ^-^^^ = {<P{x))H,t. (2.4) 

Then the generating functional r[M, t] of the one particle irreducible vertex 
functions is defined through the Legendre transformation 

T[M,t] + W[H,t] = f dxH{x)M{x), (2.5) 

with 

, , 6T\M,t] , , 



The objects of interest are the Schwinger functions (or correlation functions). 
Define 



6t{yi) . . .6t{yL)6M{xi) . . .6M{xn) 



(2.7) 



This is a functional of M and t (it is a function of its remaining arguments). 

At some value of the bare mass, the renormalized theory is massless. Set- 
ting, then, mo to this critical value, and letting t{x) become independent of 
X, t becomes a measure of the deviation away from the massless theory. If the 
source H is also independent of x, then the r*^^'^-* become functions (rather 
than functionals) of all their arguments. This is the situation henceforth 
assumed. 

Power counting yields the (primitive) degree of divergence of an individual 
graph. In less than four dimensions 

r(o,2)^r(o,4)^^^p(i,2) 

are divergent. Renormalization of these three functions gives the renormal- 
ized mass mn, the quartic coupling qr, the field strength renormalization and 
the composite field strength renormalization (Za and Z(2) respectively). The 
composite field strength is renormalized by defining {4>'^)r as Z(2) x (0^)^. 
Thus 

i<p')R = %^0^ (2.8) 

and the divergence of Z(2) renders {4>'^)r finitep|. 

In four dimensions there appears an additional divergence due to 

r(2,o)_ 

The corresponding diagram has no external legs and can never appear as 
a subdiagram. The renormalization of F^^'"^ is therefore accomplished by 
subtracting its divergent part. This subtraction does not effect the other 
p(i,A^) 'jj^g relationship between the bare and the massless renormalized 
theory is 0, |] 

^R'^\qi, ■ ■ ■ , qL;pi, ■ ■ ■ ,pn; gR, fj') = \^) z^^^ 



X {rw ^(gi, . . . , gL; Pi, • • • ,P7v; go, A) 



^2 = 1^2 



(2.9) 



in which A is the uhra- violet cutoff and fi is an arbitrary non- vanishing mass 
parameter. 

Following the renormalization prescription of 0, one can then proceed 
to expand around this critical (massless) theory. This allows one to examine 
the approach to criticality from within both the symmetric and the broken 
phases. 

The renormalization group equations express the invariance of the physics 
under a rescaling of the mass parameter fi —>■ /i(A) defined by 



yU(A) = A/i. 



(2.10) 



Following H one finds the RGE 



where 



xrSi'-"'(,ip;i(A),Af(A),jfi(A),fi(A)) = ei,A,(A) 

Sno 



Ol,nW 



-t(A)(^-^)T(to(A)) 



(2.11) 



(2.12) 



(2-L)! 

for L < 2 and zero otherwise. We will refer to this term as the inhomogeneous 
term. It comes from the need to (additively) renormalize r*^^'°^ in rf = 4. Since 
p(2,o) jg j^Q^ divergent for rf < 4, the inhomogeneous term in ( |2.11D does not 
give rise to any singular contribution there. The flow equations are 



aA,,(A) 

A dt{\) 

t(A) dX 
A dM{\) 

M{\) dX 



BigRiX)) 
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'^(gnW) 



with gn{l] 
with t(l] 



r]{gR{X)) with M(l) 



gR, 

t, 

M 



(2.13) 



and 



r{g 



R) 



.^{g- 



R) 



/^ 



d_ 

dfi 



^(2) 



.(2,0) 



rLare (?> -Q'^ £/0, A) 



q2 = |^2 ■ 



(2.14) 



The content of the RGE ( |2.11| ) is, then, that if /i is rescaled by a factor A 
then the response of qr, t and M is governed by ( p.l3| ). 

The Callan Symanzik beta function ||^ is denoted by B in ( f^.l3D . The 
form of this function can be calculated perturbatively in the renormalized 
d = (4 — e) dimensional theory. Eq. (|2.13|) then gives the behaviour of the 
running coupling constant gniX). It turns out that in order to remove the 
cutoff, the running coupling constant Qr^X) has to approach the infra-red 
(IR) fixed point g}j^ ~ 0{e). This fixed point then governs the critical region. 

The region of interest can be divided into 

(a) t > 0, H = 0. In this region there is no magnetization ( M = ). This 
is a symmetric theory. 

(b) t > and H ^ ot t < and any H. This is the region of broken 
symmetry ( in which M ^ 0). 

The RGE ( [2. Ill ) holds throughout regions (a) and (b). The point in using 
vertex functions is that it is nowhere necessary to state whether or not the 
symmetry is broken (be it spontaneously or explicitly) 0. 
The solution of ( pTl| ) is 



^{L,N) 



{q- p; t, M, gR, /i) = Z{\, ^ij(A)) 



xr(J^'^)(g;p;t(A),M(A),^^(A),^(A)) + n^,^(A), 



(2.15) 



where 



Z{X,gR{X)) 



-J 
e "R 



srW 



vi9R)+L 
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/M(A)y^ ftixy 



(2.16) 



and the inhomogeneous term is 



^d\' 



IIl,nW= —eLA^')Z{X',gR{X')). 



A 



(2.17) 



Because of the local nature of the RG, the renormalization constants of the 
infinite volume theory render finite the finite volume theory too P, ^j . We 
denote by Tj^' {q]p]t, M, gR, ji^l) the renormalized Schwinger function of 
the finite volume theory, where I denotes the linear extent of the system. The 



RGE obeyed by this Schwinger function is the same as ( |2.11| ). Its solution is 
(from (|i|)) 

^r'^\(1; P; t, M, gR, /i, /) = Z{X, gR{X)) 

xT^i'''\q;p;t{\),M{\),gni\),f^\J) + Ul,nW. (2.18) 

To prepare for dimensional analysis, we implicitly replace g^ by fi^gR to keep 
it dimensionless. Then applying dimensional analysis to the homogeneous 
term on the right hand side gives 

4^'^)(g;p;t,M,^^,/i,0 = Z(A, WA))/f (^-2)+2L-'^ 

xr^f'^) {lq;lp;lH{X),l'^M{X),gR{X),lfiX,l)+UL^N{X). (2.19) 

Since A is still at our disposal (as long as it is small enough so as to remain 
in the critical region), we choose 



Then, 



IfiX = 1. (2.20) 

^R'^\(i;p;t,M,gR,fi,i) 

xr-)(..,...(y,.^M(l),..(l),i,i) 

+I[l,n (f] . (2.21) 

In less than four dimensions, and in the critical region, the flow equations 
give [0, § 

t(A)=tA(2-^) (2.22) 

M(A)=MA-^^. (2.23) 

If Ifi is large enough then gnij-) is close to g^. Therefore, at zero momentum, 

xri,^'^) (O; 0; fi^-Hh, fiH^M, g^, 1, l) + IIl,n if] , (2.24) 



where 

/3 = ^(f/-2 + r/). (2.25) 

If fi is fixed, then 
r(f'^)(0;0;t,M,^^,/.,/) = /^+^-^F(^'^) {lh,l^M)+U,,j, (^1) , (2.26) 

where FJ^^''^^ is an unknown function of its arguments. Eq.( |2.6| ) can be 

apphed to this form for F^j' to express the external field H in terms of M. 
This gives 

4^'-)(0;0;t,i7,,^,/.,0 =/---^F;(^'-) {lhl^H)+U,,, (l) . (2.27) 



Here 6 is the usual odd critical exponent defined by 

^ d + 2-r] 



(2.28) 



d-2 + r] 
Eq. (|2.27|) is sufficient to derive the usual FSS relations in less than four 



dimensions \Wi, H, 11 



For example, the zero field susceptibility is given by 

XT\t) = fJ°'') (0;t,0,^H,/i) = f-^-'F^ilh), 

where F^ is, again, an unknown function. Putting t = then gives the FSS 
behaviour of the zero field susceptibility at the infinite volume critical point, 
Xi{t = 0)^l'-^. 

In the four dimensional version of the theory there appear certain sub- 
tleties which are not present below four dimensions. This is because the IR 
fixed point of the Callan-Symanzik function B{g^ moves to the origin as 
the dimension becomes four. Secondly, in contrast to the d < A dimensional 
case, the fixed point is now a double zero, responsible for the occurrence of 
logarithmic corrections. 

A third difference between the cases of d < 4 and d = A comes from 
the inhomogeneous term in the RGE. The graph responsible for this term 
is not in fact divergent when d < 4. Singular behaviour in less than four 
dimensions comes from the homogeneous term. In ci = 4 the inhomogeneous 



term contributes to the leading singular behaviour too. The first term re- 
mains singular however, and is responsible for divergences such as that in 
the susceptibility. 

Eg. ( p. 271) , from which the FSS behaviour of the model can be derived 
below four dimensions, was established with the help of the approximation 
gR{-^) — g*ji for large ^il. As pointed out by Brezin in p, this approximation 
fails in four dimensions. The reason is that g*^^ then becomes zero, and one 
is left with the mean field theory. 

In c? = 4, we then have to rely on a perturbative expansion in g^. To 
lowest order, the functions B{gji), rj{gR)^ ^{gn) and T{gR) are[0, |] 

Sign) = lgl-^^gi + 0{g'^), (2.29) 

v{9r) = ^9l + 0igl), (2.30) 

^ 2-lgR + Oigl), (2.31) 



r{gR) = l + 0{gn). (2.32) 

Putting yU = 1 for simplicity, these perturbative solutions, together with the 
flow equations ( |2.13|) , give for ( p.21|) 

rr'(,;P;*.M,...M)-(3^)'^— 



The coupling constant for the Schwinger function on the right hand side is 
g^ for large /. Since this is small, perturbation theory may be applied to 

calculate Fj^' 0. This gives 

rg'°)(t,M,<7^,i,0 

= CitM^ (In/)-'/' + C2M^(ln/)-' + C3t'(ln/)'/=^ (2.34) 



where ci, C2 and C3 are constants. Applying ( |2.6| ) to this yields for the external 
field 

H (t, M, gn, 1, /) ~ atMiln l)'^/^ + cgM^Qn /)-\ (2.35) 

where, again, C4 and C5 are constants. 

The free energy per unit volume in the presence of an external field is 

Wi{t, H) = MH{t, M- 1) - rg'°^(t, M; I). (2.36) 

Eqs.(g;^ and (^) give, then. 



W^Kt, H) = c;^ + 4^ + C3t^(ln/)V3, (2.37) 

where c'^ and Cg are constants and M is related to if through ( p.35| ). This 
expression is the basis of all the FSS relations derived below. 
If H vanishes, then all of the solutions of ( p.35| ) lead to 

Wi{t,Q)o^t^{\nlY . (2.38) 

3 Lattice 0^ theory and the zeroes of the par- 
tition function 

Within the path integral formulation of quantum field theory there are two 
complimentary approaches. The first is perturbation theory (in the quartic 
coupling g^. Indeed this is the basis for the considerations at the end of 
the previous section. The second approach is intrinsically non-perturbative. 
It involves the use of stochastic techniques to calculate the path integrals. 
Apart from statistical errors numerical approaches are exact, but limited to 
finite lattice volumes. 

We have used such a numerical approach — the Monte Carlo (MC) 
method — to further study the logarithmic corrections involved in four di- 
mensions. In particular, we present numerical evidence of the validity of the 
FSS formulae presented in the last section. Thus we have two independent 
approaches, whose agreement leaves little doubt that this FSS analysis in- 
deed correct. This provides support for the validity of the analyses presented 
in , and [Q , and for the triviality of 04 theory. 

10 



The usual regularization for a numerical approach replaces the space-time 
continuum by a lattice. This is, of course, entirely equivalent to the use of 
the momentum cut-off in sect. 2. We use a regular hypercubic lattice of unit 
intersite spacing. If t is independent of x in ( p.l|) , the lattice parameterized 
action in the absence of a source field and with finite differences replacing 
derivatives reads 

-«5:</..<^.+^+^0^^A^ (0^.-1)' . (3.1) 

X,fJ, X X 

Here the hopping parameter k and the quartic coefficient A correspond, in a 
sense, to the mass and quartic coupling of the continuum theory respectively. 
Taking A to infinity gives the Ising limit of the model. Here, the fields (px 
take only values from the set {±1}. The universality hypothesis, which comes 
from experience in statistical physics, implies that no information should be 
lost in going to the Ising extreme. I.e., the Ising model and the 0^ model 
with arbitrary k and A should be in the same universality class and exhibit 
the same scaling behaviour (for a related MCRG study cf. ||12|). The vacuum 
to vacuum transition amplitude of the quantum field theory becomes the 
partition function of the Ising model. 

In the presence of an external field the Ising model can be defined by the 
partition function 

Zin,H) = ^j:e^'^''' (3-2) 

where 

-5 = E E <^-'/'-+M , M = E0x . (3.3) 

X /x=l X 

Here, the Boltzmann factor has been absorbed into the hopping parameter 
K, and into the reduced external field h = k x H. The sum runs over all J\f 
possible configurations of the spin field on the d dimensional lattice, and the 
normalization ensures Z(0,0) = 1. We may reexpress the partition function 

by 

N dN 

Z{K,h) = Y. Y. p(^,M)e"^+^^ 

M=-N S=~dN 

N dN 

= E p(«:;M)e'*^= E PiS;h)e'^' , (3.4) 

M=-N S=-dN 

11 



where A^ is the number of sites on the lattice. The spectral density p{S, M) 
denotes the relative weight of configurations having given values of S and M. 
By p{k; M) and p{S] h) we denote the correspondingly integrated densities. 

Z is a polynomial in the fugacity e^^ (degree A^) and e^'^ (degree dN/2). 
The coefficients of the polynomial in e^^ for real constant k are real and pos- 
itive, as are those of the polynomial in e^'^ for real constant h. A knowledge 
of the zeroes of the partition function is equivalent to a knowledge of Z itself 
(and of all functions derivable from it). In particular, the critical behaviour 
of Ising-type systems can be analysed through its partition function zeroes 
instead of more traditional methods involving real parameters. 

The study of partition function zeroes in general was initiated by Yang 
and Lee in 1952 |]T3. The Lee- Yang theorem states that for ferromagnetic 



systems all of the zeroes of the partition function in the external ordering 
magnetic field variable lie on the imaginary axis for real temperatures. Fisher 
was the first to analyse the zeroes in the complex temperature (or mass) plane 



Ij]. Thus we refer to partition function zeroes in the temperature plane as 



Fisher zeroes and to those in the complex plane of external fields as Lee- 



Yang zeroes. With the exception of systems which are self dual []15[, there 
exist no simple general results of the Lee- Yang type concerning the locus of 
Fisher zeroes. Thus the vast majority of studies have been of a numerical 
nature (see, however |T^, |T^ and references therein). 



Itzykson, Pearson and Zuber [T^ initiated the study of FSS of partition 



function zeroes. Their analysis was confined to less than four dimensions 
with power-law scaling behaviour and corrections. This was later extended 



to dimensions above (not including) four in [0. The latter is also restricted 
to purely power-law scaling behaviour. In this section, the corresponding 
FSS theory is presented for four dimensions where logarithmic corrections 
are manifest. 

Denote by Ci(t) and Xiif) the specific heat and magnetic susceptibility 
(per unit volume) respectively of a system of linear extent / and at a reduced 
temperature t in zero external magnetic field (cf. sect. 3). Twice differentiat- 
ing the free energy in the perturbative RG formula ( [^.38D gives 

G(t) oc (In/)^. (3.5) 

The total free energy at the critical temperature in four dimensions in 



12 



the presence of an external field is given by ( p.37| ) as 

l\\nl)-^Hl. (3.6) 

The partition function is therefore 

Zi{t = 0,H) = Q{l\\niy^Hl). (3.7) 

If at some (complex) value of H the partition function vanishes, then, for 
this value of H, 

hI =l-\\nl)-^^Q-\0). (3.8) 

Therefore 

Hj ocl-%\nl)--i (3.9) 

where the constant of proportionality depends on the index j of the zero. 
This is the FSS formula for Lee- Yang zeroes in four dimensions. 
Eq.( p.38|) gives for the total free energy (when H = 0) 



Fi(t,if = 0) oc/V(ln/)3. (3.10) 

The partition function is the exponential of this, i.e., 

Zi{t, H = 0) = R (iH^ (Int)'/') . (3.11) 

If R vanishes, then, 

t^l^ {Intf^ = Rj\0) (3.12) 

where j indicates the index of the zero. Therefore, for the j^^ zero, 

t,oc/-2(lnO"'/' (3.13) 

where the proportionality constant depends, again, on j. 

The scaling relations for the partition function zeroes can be used to 
find the behaviour of the thermodynamic functions as well. The partition 
function is a polynomial and as such can be written in terms of its zeroes. 
Let Hj be the j^^ Lee- Yang zero for a system of linear extent l. Then, the 
partition function is 

Zi{k,H)ocI[{H-Hj). (3.14) 

j 
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The magnetic susceptibility is given by the second derivative of the Gibbs 
free energy with respect to H. This gives, in d = 4, 

Xi{i^, H)(x-Y. 777 7772 • (3-15) 

Therefore the susceptibihty at the critical value of H (namely at H = 0) is 

1 1 

X^/^^O) oc -^— 2- (3.16) 

* j ^J 

Eq.( |3.9| ) then gives the FSS formula for the zero field susceptibility in four 
dimensions as 

XiiK,,H = 0)^l^\nl)K (3.17) 

A similar calculation for the Fisher zeroes leads to the recovery of the 
FSS formula for specific heat ( |3.5| ). Let Kj be the j^^ Fisher zero for a system 
of linear extent /. Then, in zero field, the partition function is 

Zi{k)ocY[{k-Kj). (3.18) 



The specific heat is given by the second derivative of the free energy with 
respect to t. This gives, in d = 4, 

C/(«:) = -^E-^^- (3-19) 

Therefore the specific heat at the critical value of k is 

where Tj is the 'reduced' position of the j^^ zero: 

Tj —— t\jj ^c* (O.^-L) 

Eq. (|3:T3|) gives 



r. 



-J oc/-2(ln/)"6 (3.22) 



14 



where the constant of proportionahty depends on the index j. Thus the FSS 



formula for the specific heat in four dimensions is [20| 



1 



Q(kJ oc (ln/)3 . (3.23) 

The partition function zeroes provide an alternative way to view the onset 
of criticality. As the system size increases towards infinity, the zeroes tend to 
pinch the real H or k, axes ( ( |3.9|) and ( p. 13 )). Thermodynamic observables 



such as the specific heat and magnetic susceptibility become divergent. This 
applies to the correlation length as well. The FSS formula for the correlation 
length of a four dimensional system also involves logarithmic corrections. 
This was derived by Brezin |^ for a system of extent / in all directions. At 
the infinite volume critical point k, = Kc, one has 

e/(fi:e)oc/(ln/)i (3.24) 

This suggests that a FSS variable should indeed be defined by 

^oo(/t) t'2 \ Int I 6 



eK«^c) l{lniy 



(3.25) 



in four dimensions P]. 



4 Numerical calculations and results 

We now want to report on our numerical calculations which confirm the scal- 
ing picture of sect. 2. In particular, we want to identify the multiplicative 
logarithmic corrections to FSS. Such logarithmic corrections have been no- 
toriously difficult to verify numerically (see e.g. M and pOl). However, the 



advent of more efficient cluster algorithms ||21|, |2^ has greatly improved the 
quality of Monte Carlo calculations for bosonic spin systems like the Ising 
model. We suggest — and the quality of our results supports our proposal — 
that a study of the FSS of partition function zeroes lends itself more readily 
to the detection of logarithms than do the more traditional thermodynamic 
quantities such as specific heat. 

The first numerical calculations of partition function zeroes appeared in 



the 1960's ||2^. Such early work involved exact calculations of the density of 

15 



states (spectral density) p{S, M) and were therefore confined to very small 
lattices. (See [^ for a list of references and early history). 

The next major step concerning numerical calculations was made by Fal- 
cioni et al.[^ and by Marinari et al.[Q in the early 1980's. They were the 



first to use approximations to the density of states in the form of histograms 
to study critical phenomena. It is clear that straightforward analytical con- 
tinuation can take ( p.4|) to the complex k or H plane. Thus the histogram 
technique can be used for a precise numerical determination of the complex 
partition function zeroes. 

Numerical methods received a further boost with the development of 
techniques whereby a number of Monte Carlo constructed histograms can 
be combined to form one 'multihistogram' [E^, ESI. These provide a better 



approximation to the spectral density over a wider range of the parameter k 
or H. 

We now present some details of our own numerical calculations. The data 
were taken on lattices of size from 8^ to 24^ using the Swendsen-Wang cluster 
algorithm. Histograms were determines at /i = and at various values of 
K, close to the pseudocritical one (chosen to be that value of k where the 
specific heat peaks). Table 1 provides a list of lattices sizes, the values of k 
at which the simulations took place, as well as a summary of the statistics. 

For the determination of p{S; h = 0) the various 'raw' histograms were 
suitably combined following [^]; no binning was used. With (|3^ ) this allows 
one to construct Z{k, h = 0) in the complex neighbourhood of the real k- 
values and to determine nearby Fisher zeroes. 

For p{k,; M) we binned each of the raw {S, M)-histograms in a 256 x 256 
array and then combined for each M-bin the corresponding S'-subhistograms 



according [^. This then allows us to obtain an optimal p{ko; M) for arbi- 
trary K,Q in the considered domain. From this Z{ko, h) may be determined for 
not too large values of (imaginary) h. Below we present results for Lee- Yang 
zeroes evaluated at kq = n^- As a consistency check we also determined the 
zeroes coming from a single M— histogram corresponding to a simulation at 

The errors in the quantities calculated from the multihistograms were 
estimated by the jackknife method, i.e. the data for each lattice size were 
cut to produce 10 subsamples leading to different multihistograms and thus 



to different results, whence the variance and bias were calculated f2£ 



From the Lee- Yang theorem |jT3| it is known that the zeroes in H all lie on 
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the imaginary axis for any lattice size /. The search for the Lee- Yang zeroes 
is therefore technically easier than for the Fisher zeroes. In the later case we 
used a Newton-Raphson type algorithm. In order to avoid instabihties due 
to the large numbers involved, and since Zi{k) never vanishes for real k, the 
Fisher zeroes were in fact found as local steep minima in \Zi{K)/Zi{Iie k)\ . 

We now come to the FSS analysis of the Fisher zeroes. The positions of 
the closest two Fisher zeroes obtained from the multihistograms are listed in 
table 2 (where Kj represents the j*^ Fisher zero). Since we can confine the 
scaling analysis to the imaginary parts of the zeroes we avoid the necessity of 
knowing the infinite volume critical value of k. In fig. la we plot the logarithm 
of the imaginary part of the position of the first Fisher zero against the 
logarithm of the lattice size /. A linear fit to the slope (— -) gives u = 0.479(1) 
which is slightly below the mean field value of |. This deviation from the 
mean field value is due to the presence of logarithmic corrections, which we 
have neglected in this first fit. A corresponding analysis applied to the second 
Fisher zeroes gives u = 0.467(8). 

Assuming that the leading scaling behaviour is indeed proportional to 
/^, we can proceed to search for multiplicative logarithmic corrections. To 
this end, we plot in fig. lb ln(/^ImKi) versus In (In/). A negative slope is 
clearly identified and is in good agreement with the scaling prediction of 
— |. In fact, a fit to all five points gives a slope —0.217(12). Excluding the 
point corresponding to / = 8 gives a slope of —0.21(4). The solid line is the 
best fit to the points corresponding to / = 12 ... 24 assuming the theoretical 
prediction — | from ( p.l3|) . 

The errors in the second (and higher index) Fisher zeroes are too large 
to warrant a corresponding analysis. 

Now that the logarithmic corrections to the FSS behaviour of the Fisher 
zeroes have been established, we may proceed to determine the infinite vol- 
ume critical hopping parameter Kc from 



,-2- ' 



\kj - Kc\ (xl (In/) 6 . (4.1) 

Using the first Fisher zeroes, we find Hc — 0.149703(15) in good agreement 



with the value 0.149668(30) from high temperature expansions [^ . 

Mean field theory [|I^] predicts that the angle ip at which the Fisher zeroes 
depart from the real axis should be ^. There exists, unfortunately, no FSS 
theory for this quantity. We list, in table 3, our measurements of this angle 
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defined by the first and second Fisher zeroes and the real axis for the five 
lattice sizes analysed and plot it in fig. 2. The average value compares well 
with the mean field prediction. 

Let us now discuss the results for the Lee- Yang zeroes. For all real k 
they have to lie on the imaginary /i-axis. At Hc they should scale according 
to the FSS formula ( p.9| ). Table 4 lists the positions of the first two Lee- Yang 
zeroes as obtained from the multihistogram at Kc = 0.149703. Fig. 3a is a log- 
log plot of the imaginary parts of the positions of the first Lee- Yang zeroes 
against the lattice size. The resulting slope is —3.083(4); this compares well 
with mean field prediction of —3. The deviation, again, may be explained 
by logarithmic corrections. 

To identify the logarithmic corrections, we plot in fig. 3b In {iHmhi) against 
In (In /). For the first Lee- Yang zeroes a best fit to all five points gives a slope 
of —0.204(9) which compares well with the theoretical prediction of —^ from 
( p.9| ). Excluding the smallest lattice, a best fit to the remaining four points 
gives a slope —0.22(3). The solid line in fig. 3b is the best fit to the last four 
point with given slope —^. 

The errors in the positions of the second nearest Lee- Yang zeroes are 
only about twice that of the corresponding first index zeroes and these can 
also be used to analyse FSS. The results again are good agreement with the 
expected scaling behaviour. 

5 Conclusions 

We have used RG techniques to derive the FSS behaviour of the cp'^ theory in 
d=4, placing particular emphasis on the partition function zeroes. These for- 
mulae were then tested in a non-perturbative fashion — with high precision 
numerical methods. Of primary interest are the multiplicative logarithmic 
corrections to the leading power law scaling behaviour. These logarithmic 
corrections are clearly identified from the scaling behaviour of the closest 
Fisher and Lee- Yang zeroes and are in good quantitative agreement with 
the theory. Higher index partition function zeroes and thermodynamic ob- 
servables such as the specific heat exhibit logarithmic corrections too. 

The high precision of the numerical results and the good quantitative 
agreement with the analytical predictions of sect. 3 provide non-perturbative 
evidence for the existence of a double zero in the Callan-Symanzik beta 



function. The fixed point responsible for non-trivial behaviour in d < 4 di- 
mensions has moved to the the position of the trivial fixed point in d=4. The 
results support the assertion that 0^ theory is trivial in four dimensions. 
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Tables 

Table 1: For each lattice size L we list the values of k at which MC simu- 
lations were performed. The corresponding h value is always zero. In paren- 
theses we give the number of measurements in units of 1000 (^^/lOOO). 



L 



Ko(# /lOOO) 



0.149709 (400), 
0.1525 (400), 
0.1533 (400), 



0.1506 (200), 
0.1527 (400), 
0.1540 (200) 



0.1515 (200), 0.1520 (200), 
0.1529 (400), 0.1531 (400), 



12 



0.149709 (200), 
0.1510 (200), 



0.1498 (100), 0.1503 (100), 0.1508 (200), 
0.1512 (200), 0.1514 (200), 0.1520 (100) 



16 



0.1492 (50), 
0.1503 (100), 
0.1506 (150), 
0.1518 (50) 



0.1496 (50), 0.149709 (150), 0.1500 (50), 
0.1504 (100), 0.1505 (150), 0.15054 (150), 
0.1509 (100), 0.1511 (50), 0.1513 (50), 



20 



0.1495 (30), 0.149709 (80), 0.1498 (30), 0.1499 (50), 

0.1500 (50), 0.1501 (100), 0.1502 (80), 0.1503 (80), 

0.1504 (50), 0.1507 (30) 



24 



0.1495 (20), 
0.1500 (80), 



0.149709 (80), 
0.1501 (80), 



0.1498 (108) 
0.1502 (80), 



0.1499 (68), 
0.1504 (20) 



Table 2: The positions of the first and second Fisher zeroes as obtained 
from the jackknifed multi-histograms. 



L 


Re(/ti) 


lm(fi:i) 


Re(/t2) 


lm(/t2) 


8 


0.152156(10) 


0.004046(10) 


0.154195(104) 


0.006085(127) 


12 


0.150802(7) 


0.001733(10) 


0.151652(86) 


0.002615(137) 


16 


0.150322(7) 


0.000948(5) 


0.150913(38) 


0.001352(41) 


20 


0.150095(5) 


0.000595(7) 


0.150397(55) 


0.000875(36) 


24 


0.149972(3) 


0.000414(5) 


0.150198(40) 


0.000574(48) 
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Table 3: The angle ipi^2 between the first and second Fisher zeroes as ob- 
tained from the jackknifed multi-histograms. Mean field theory predicts 

(^ = 7r/4 ~ 0.785. 



L 


V'1,2 


8 
12 
16 
20 

24 


0.785(62) 

0.804(139) 

0.599(89) 

0.747(175) 

0.617(246) 


mean 


0.731(45) 



Table 4: The positions of the first two Lee- Yang zeroes as obtained from 
the jackknifed multihistograms at /t = 0.149703. The real part of the zeroes 
is always zero. 



L 


lm(/ii) 


Im(/i2) 


8 
12 
16 
20 

24 


0.0022294(32) 

0.0006367(23) 

0.0002637(9) 

0.0001327(7) 

0.0000749(5) 


0.0049488(71) 
0.0014111(22) 
0.0005845(33) 
0.0002949(18) 
0.0001656(7) 
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Figures 

Fig. 1: (a) The imaginary part of the Fisher zeroes closest to the real /t axis 
vs. the logarithm of the lattice size L. The straight line is a fit to all points 
with slope -2.088(6), slightly differing from the expected value —1/v = —2. 
This difference is due to logaritmic corrections to scaling as shown in (b), 
where we plot ln{L'^lm.Ki) vs. InlnL. 

Fig. 2: The impact angle ipi^2 of the closest two Fisher zeroes, compatible 
with the expected value 7r/2 (full line) within the error (shaded area). 

Fig. 3: (a) The imaginary part of the Lee- Yang zeroes closest to h = 
vs. the logarithm of the lattice size, for k = 0.149703. The straight line is 
a fit to all points with slope -3.083(4), slightly differing from the expected 
value —1/6 = —3. This difference is due to logaritmic corrections to scaling 
as shown in (b), where we plot ln(L^Im/ii) vs. InlnL. 
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